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Abstract. In this paper, a new general decomposition theory inspired from modular graph decompo- 
sition is presented. Our main result shows that, within this general theory, most of the nice algorithmic 
tools developed for modular decomposition are still efficient. 

This theory not only unifies the usual modular decomposition generalisations such as modular de- 
composition of directed graphs or decomposition of 2-structures, but also star cutsets and bimodular 
decomposition. Our general framework provides a decomposition algorithm which improves the best 
known algorithms for bimodular decomposition. 

1 Introduction 

Several combinatorial algorithms are based on partition refinement techniques [20]. Graph algo- 
rithms make an intensive use of vertex splitting, the action of partitioning classes between neigh- 
bours and non-neighbours of a vertex. For instance, all known linear-time modular decomposition 
algorithms [14, 18, 8, 16, 6, 4, 13, 9] use this technique. 

In bioinformatics also, the distinction of a set by a element, called a splitter, seems to play an 
important role, as for example in the nice algorithm of [21], which computes the set of common 
intervals of two permutations. 

In this paper we investigate an abstract notion of splitters and propose a dual formalism based 
on the concept of homogeneity. Our aim is a better understanding of the existing modular decom- 
position algorithms by characterising the algebraic properties on which they are based. Our main 
result is that, within this general theory, most of the nice algorithmic tools developed for modular 
decomposition [14, 18, 8, 16, 6, 4, 13, 9] are still efficient. 

This theory not only unifies the usual modular decomposition generalisations such as modular 
decomposition of directed graphs [17] or decomposition of 2-structures [10], but also allows to handle 
star cutsets, and the bimodular decomposition [11]. Notice that our general framework provides a 
decomposition algorithm which improves the best known algorithm for bimodular decomposition. 

The paper is structured as follows: first we detail this new combinatorial decomposition the- 
ory, then we present a general algorithmic framework, and we finish by listing some interesting 
applications. 



2 Homogeneity, a new viewpoint 

Throughout this section F is a finite set. The family of all subsets of V is denoted J'(l^). An 
reflectless triple is {x,y,z) C with x y and x ^ z. Reflectless triples will be denoted {x\yz) 
instead of (x, y, z) since the first element does not play the same role. Let H he a relation over the 
reflectless triples of V. Given s G y, the relation is a binary relation on V defined as Hs{x,y) 
if and only if H{s\xy). 



Definition 1 (Homogeneous relation). H is a homogeneous relation on V if, for all s E V, 

Hg is an equivalence relation onV\ {s}: 

- (Symmetry): y s,x,y & V, H{s\xy) <^ H{s\yx). 

- (Refiexivity): y s,x E V, s^x =^ H{s\xx). 

- (Transitivity): y s,x,y, z E V, H{s\xy) A H{s\yz) H{s\xz). 

Definition 2 (Homogeneous sets). Let H be a homogeneous relation. X CV is homogeneous 

with respect to an clement s ^ X if H{s\xy) for all x,y £ X. If X is not homogeneous w.r.t. s 
then s distinguishes X, or is a splitter of X. Let §x be the set of all elements distinguishing X 
and s{X) = |Sx|- 

M C.V is a homogeneous set if M ^ ^ and for all x not in M , M is homogeneous w.r.t. x. 
In other words s{M) = 0. The family of homogeneous sets for a homogeneous relation H on V is 
denoted 3^h or 5" if not ambiguous. 

Remark 1. From the definition it is obvious that, given a homogeneous set M, if -'H(^s\xy) for some 
x,y e M then s G M. 

Two sets A and B overlap 'd An B, A \ B and B \ A are all nonempty. It is denoted A^B. 
The symmetric difference of two sets A and B, denoted AAB, is {A\ B) U {B \ A). Let us now 
enumerate some properties of the homogeneous relations and sets. 

Proposition 1. If s distinguishes X then ^H{s\xy) for some x,y £ X. 

Thanks to this proposition, it is exactly equivalent to define the homogeneity relations as ternary 
relations on reflectless triples - as presented here - or as relations between elements and subsets 
of V. A homogeneous relation H{s\X) for X C F and s \ X simply is a transitive relation: if 
A(^B and H{s\A) and II{s\B) then H{s\A U B). The relation H has the same properties than H 
and the homogeneous sets are the same. We find the ternary relation much simpler. 

Proposition 2. For all A, B e 3^ if A<^ B then (A n S) G and (A U 5) G 5'. 

This property is called closure under intersection and union. It is easy to check and can be used 
to prove: 

Proposition 3 (Lattice structure). Let H he a homogeneous relation on V and 3^'^ = J'jy U{0}. 
is a lattice. 

Proof. Since G and thanks to Proposition 2, the intersection of any two elements of 
belong to 3"^. It is the infimum of two sets, since any member of contained in both A and B 
is contained in AD B. Let us consider the family S of all set of containing both A and B. It 
is nonempty (V is a member). Since 3"^ is closed under intersection, S admits a unique smallest 
(w.r.t. inclusion) element, the intersection of all its members, that is the supremum of A and B. □ 

This lattice is a sublatticc of the boolean lattice (hypercube) on V. Moreover, if we consider 
Ae3 such that |^| > 1, 3iA) = {3 e 3h and 3 ^ A} then {3{A), C) is a distributive lattice. Let 
us now define some useful types of homogeneous relations. 

Definition 3. A homogeneous relation H fulfills 
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- Al if\/x,y,zeV,H{x\yz) A H{y\xz) =^ H{z\xy). 

- A2 ify s,t,x,y eV, H{x\st) A H{y\st) A H{t\xy) H{s\xy). 

- A3 if\/ s,t,x,y €V, Hlx\st) A Hiylst) A Hlt\sx) A H(t\sy) H(s\xy). 

- A4 ifWx,y,zeV,-'H{x\yz) A -'H{y\xz) ^ H{z\xy). 

Proposition 4 (Quotient). Let H be a homogeneous relation. Then, H satisfies A2 if and only 
if for all homogeneous set M , for all x,y G M and s,t G V \ M , H{x\st) H{y\st). 

It is a simple rewriting of A2 but it enlightens that, from a homogeneous set, one can pick a 

representative element. Indeed, elements in a homogeneous set M uniformly perceive a set X 
not intersecting M: if one element of M distinguishes X then so do all. This allows to shrink a 
homogeneous set M into a single element, the quotient by M. 

Given X CV one can define the induced relation H[X] as H restricted to reflectless triples of 
X^. If X is a homogeneous set we have the following nice property: 

Proposition 5 (Restriction). Let H be a homogeneous relation, M a homogeneous set and M' C 
M. M' e S'him] ^ M'eS^H. 

Recursivity can therefore be used when dealing with homogeneous sets. Notice that the propo- 
sition is not always true if M is not a homogeneous set. The Quotient and Restriction properties 
were used first with modular decomposition and are useful for algorithmic [19]. 

3 Submodularity of homogeneous relations 

Definition 4. A set function ^ : ^ M is submodular if and only if for all X,Y C V 

fi{X) + /i(y) > n{X UY) + n Y) (see e.g. [12]). 

Theorem 1 (Submodularity). Let H he a homogeneous relation. The function s counting the 

number of splitters (defined, in Definition 2) is submodular. 

Proof. Let us for convenience define s(0) = —\V\. It suffices to prove s{X) + s{Y) > s{X U y) + 
s{X DY) for all overlapping X,Y C V. So let X,Y C V he two overlapping sets. Obviously, 
S>xnY = (Sxny \ Y, Sxnr H Y) . 

As an element distinguishing X does not belong to X, the partition §xur = (Sxuy \Sx; Sxur H 
§x ) can be reduced to §xuY = (Sxuy \ Sx , SixXiX UY)). Similarly, Sy = (Sy \ §xnY , Sxny \ . 
Finally, §x = (Sx \ Y, (Sx nY)\ §xnY, (Sx n F) n §xny) can be reduced to $x = (Sx \ U 
Y),{§xnY)\ Sxny , Sxny n Y). Hence, |§x| + |Sy | - |§xuy I - |§xny I = |(§x n F) \ §xny I + |Sy \ 
§xny I — |Sxuy \ Sx|- 

To achieve proving the theorem, we prove that Sxuy \Sx ^ Sy \Sxny • Indeed, let z G Sxuy \Sx- 
Then, z ^ X UY and for all x,y e X, we have H{z\xy). 

Now, suppose that z ^ §y . Since z is not in XUY, we have H{z\xy) for all x,y eY. Furthermore, 
as X and Y overlap and thanks to the transitivity of H, we have z ^ §xuY, which is a contradiction. 

Finally, supposing z G §xnY would imply z € 8x- □ 

Submodular functions are combinatorial objects with powerful potential (see e.g. [12]). The- 
orem 1 enables the application of this theory to homogeneous relations. In [21], T. Uno and M. 
Yagiura gave a (restricted) version of this theorem, and constructed a very nice algorithm com- 
puting the common intervals of a set of permutations. It would be interesting to generalise this 
approach to any homogeneous relation, as done in [3] for modular decomposition. 
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4 Strong homogeneous sets and Primality 

In a family J" of subsets of V, a subset is strong if it overlaps no other subset of 9^. The other subsets 
are weak. Let us suppose the family 3^ contains V and the singletons {v} for every element. Then V 
and {v}v^v form the trivial strong subsets. The set inclusion orders the strong subsets into a tree. 
This is a quick proof that there are at most 2\V\ — 1 strong subsets (and at most |y| — 2 nontrivial 
ones), as the tree has no internal node of degree 1. 

The parent of a (possibly weak) subset M is the smallest strong subset Mp properly containing 
M, and M is said to be a child of Mp. If M is strong, Mp is by definition its parent in the inclusion 
tree. 

An overlap class is an equivalence class of the transitive closure of the overlap relation QD on 
3^. The support of an overlap class C = {Ci, ...Ck} is Ci U ... U Cfe. A is an atom of the overlap class 
if it is included in at least one subset Q, and it does not overlap any subset of the class, and is 
maximal for these properties. All the atoms of a class form a partition of its support, the coarsest 
partition compatible with the class. An overlap class is trivial if it contains only one subset; it is 
then clearly a strong one. 

A strong subset is prime if all its children are strong, and decomposable otherwise. It is a classical 
result of set theory that 

Lemma 1. // H is a family closed under union of overlapping sets, then there is an one-to-one 
correspondence between the nontrivial overlap classes of 3 and the decomposable strong subsets of 
3. More precisely, the overlap class C associated with a decomposable subset D is simply the set of 
weak children of D, and the support of C is D. 

The overlap class associated with a decomposable node is simply the set of weak children of 
this node. Of course we apply all these notions on homogeneous set families. On partitive families, 
the strong homogeneous sets plays a very important role since they are exactly a coding, in 0(|y|) 
space, of the possibly 21^' subsets of the family (see the upcoming section). 

Theorem 2. Let H be a homogeneous relation and Z be the family of homogeneous sets containing 
X but not y, and maximal for this property, for all x and y. The strong homogeneous sets of H are 
exactly the supports and atoms of all overlap classes of Z 

Proof. First, remark that, thanks to the closure under union of overlapping sets, the supports and 
atoms of every overlap class of Z are strong homogeneous sets. Lemma 1 tells they can not be 
overlapped by an element of Z and if one, A is overlapped by an homogeneous set B ^ Z then for 
X € A\B, the maximal homogeneous set containing y but not x overlaps A, a contradiction. So the 
family of supports and atoms is included in the family of strong homogeneous sets. Conversely, let 
us prove that if M is a strong homogeneous set then it is the support or an atom of some overlap 
class. We shall distinguish four cases. Let Mp be the strong parent of M (for M ^V). 

1. M is trivial {V or {v}). There is no problem. 

2. M is decomposable. It has k strong children Mi..Mfc. Let us pick an element Xi in each Mj. 
Then for all i and j we consider the maximal homogeneous set containing Xi but not Xj. They 
form an overlap class of Z. Its support is M, thanks to Lemma 1 

3. M is prime and Mp is prime. Then for all x G M and all y G Mp \ M M is the maximal 
homogeneous set containing x but not y. As it is strong, it belongs to a trivial overlap class and 
is equals to its support. 
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4. M is prime and Mp is decomposable. Then for all x G Mp \ M then M is included in some 
maximal homogeneous set not containing x (the one that contains the vertices of M). Let 
us consider the intersection / of all subsets of {M^ \ x € Mp \ M}. It is an atom of the overlap 
class associated with Mp and thus is strong. As M is a children of Mp, I = M. □ 

This theorem leads to a 0(|l^p)-time algorithm in Section 6.4. 
5 Partitive families of homogeneous sets 

A generalisation of modular decomposition, known from [5], less general than homogeneous relations 
but more powerful, is the partitives families. 

Definition 5. A family 3^ C J'(F) is weakly partitive if it contains V and the singletons {v} for 
all V E. V, and is closed under union, intersection and difference of overlapping subsets, i. e. 

AeS^H A Be^H A A'^B ^ AnBeJ^n A AuBe'JH A A\Be'JH 

Furthermore a weakly partitive family 3^ is partitive if it is also closed under symmetric difference: 

AeS'H A BeS^H A A(^B AABeS^H 

As mentionned before, strong subsets of a weakly partitive family 3^ can be ordered by inclusion 
to a tree. Let us define three types of strong subsets, i.e. three types of nodes of the tree: 

— prime nodes who have no weak children, 

— degenerate nodes such that all union of strong children of the node belongs to JF, 

— linear nodes such that there is an ordering of the strong children such that a union of children 
belongs to if and only if they are consecutive in this ordering. 

Theorem 3 {[5]). In a partitive family, there exists only prime and degenerate nodes. In a weakly 
partitive family, there exists only prime and degenerate and linear nodes. 

The strong subsets are therefore an 0(|F|) space encoding of the family: it is enough to type the 
nodes into complete, linear or prime, and to order the children of the linear nodes. All weak subsets 
can be output just by making simple combinations of the strong children of decomposable (complete 
or linear) nodes. Now, the following properties state that the homogeneous relations are a proper 
generalisation of (weakly) partitives families. 

Proposition 6. Let H be a homogeneous relation. If H fulfills Al , or if H fulfills A2, then H 
fulfills A3. 

Remark 2. This proposition allows to classify the homogeneous relations. There exists homogeneous 
relations fulfilling Ai but not Aj for all i and j excepted the two implications of Proposition 6. 

Proposition 7. // a homogeneous relation H fulfills A2 or A3, then 3^h is a weakly partitive 
family. 

If H fulfills A 1 , then 3^h is a partitive family. 
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Algorithm 1: SuimIIcsI Ilomogcucous set coulaiuiug S 

Let X be an element of S, M := {x} and F := S\ {x} 
while F is not empty do 

pick an element y in F ; F := F \ {y} ; M := M U {y} 

for every element z do 
|_ if H{z\x,y) then F := FLI{z} 

output M (now equals to SHS{S)) 



Proof. Let us suppose A G 3^h and B G S^h and A'^B. Thanks to transitivity an element not 
in A LI B can not distinguish AU B (it would distinguish ^4 or B). As an element not in A can 

not distinguish A and an clement not in B can not distinguish B, then no element can distinguish 
An B. For the same reason, only an element of ^ n -B can distinguish A\B oi AAB. 

If s & An B distinguishes A \ B, then this set contains x and y such that ^H{z\xy). But as 
B \ ^4 is nonempty it contains t and we have H{x\st) and H{y\st) and H{t\sx) and H{t\sy) and 
H{t\xy). Then both A2 and A3 are violated. 

Let us suppose Al holds. As A3 also holds A\B and B\A are homogeneous sets. If z G AOB 
distinguishes AAB, then there exists x & A and y & B such that -^H{z\xy). But since H(x\yz) and 
H{y\xz) -4i is contradicted. □ 

6 Homogeneous Set Algorithms 

In the following, we consider a fixed ground set V and a homogeneous relation H on V, that are 

the input of all algorithms described here. The input H consists in |y| partitions (the equivalence 
classes of for each x) and thus can be stored in 0(|Fp) space, instead of the naive 0(|Fp) 
space representation storing all triples. 

6.1 Smallest Homogeneous Set containing a subset 

Let H he a homogeneous relation on V and S be a nonempty subset of V. As is closed under 
intersection, there is an unique smallest homogeneous set containing S, the intersection of all 
homogeneous set containing S, denoted henceforth SHS{S). 

Theorem 4. Algorithm 1 computes SHS{S) in 0{\V\.\SHS{S)\) = 0{\V\^) time. 

Proof. Time complexity is obvious as the while loop runs |M| — 1 times and the for loop \V\ times. 
The algorithm maintains the invariant that every splitter of M is in F. When M is replaced by 
M U {y}, every element that distinguishes M U {y} distinguishes x from y, or already is in F. The 
algorithm ends therefore on a homogeneous set that contains S, and thus we have SHS{S) C M. 
IfM^ SHS{S) let V be the first element of M \ SHS{S) added to F (eventually added to M). It 
distinguished two elements x and y from SHS{S), contradicting its homogeneity. So SHS{S) = M. 

□ 

6.2 Maximal Homogeneous Sets not containing an element 

Proposition 8. Let H be a homogeneous relation on V and x G V an element. As is closed 
under union of intersecting subsets, there is an unique partition of V \ {x} into Si...Sk such that 
every Si is a homogeneous set of 3^h and is maximal w.r.t. inclusion in In- 
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We call MHS{x) C J'(F) this partition of Maximal Homogeneous Sets not containing x. We 
propose a partition refining algorithm [20] . It is obvious that 

Lemma 2. Every homogeneous set (especially the maximal ones) not contains x is included in a 
X- class HI. of H. 

Therefore our algorithm starts with the partition P = {H^..!!^} of the x-classes of H. Then 
the partition is refined (classes are splitted) using the following rule. Let y be an element, called 

the pivot, and C{y) the class of Partition P containing y. 

Rule. For a given pivot y, split every class of P, excepted C{y), into C fl Hy,...,C n Hy 
Notice that a class is actually splitted in many new classes iff it is distinguished by y. 

Lemma 3. Starting from the partition Pq = {H]^..H^}, the application of the refining rule (for 
any pivot in any order) until no class can be actually splitted, produces MHS{x) 

Proof. The refining process ends when no pivot can split a class, i.e when every partition class is a 
homogeneous set. Let us suppose one of these homogeneous sets M is not maximal w.r.t. inclusion: 
it is included in a homogeneous set M', itself included in a x-class iJ* . Let us consider the pivot y 
that first broke M'. It can not be out of M', as M' is homogeneous, nor within M', as a pivot does 
not break its own class. But M' was broken, contradiction. □ 

Let us now implement this lemma into an efficient algorithm. The idea for saving time is that, 
after a class is splitted by y, it has not to be examined one more time but, if the former class C 
containing y is later splitted into new classes Ci..Ca- W.l.o.g, suppose y € Ci. Then only the new 
classes C2--Ca must be examined when y is the pivot again. Every element is thus examined once 
for every pivot, leading to an 0(|yp) time complexity. This is implemented in Algorithm 2. 

We implement this idea using groups. The partition is refined from Pq into Pi then P2 and so 
on. The group of an element of Pi is its class in If the classes are implemented using a linked 

list, and if the classes are only splitted into new classes that follow consecutively in the linked list, 
then the group boundaries are simply markers in the linked list. 

A classical trick of partition refining [20, 14] is using a refining set R . Every class C of P 
can be splitted into C OR and C\R in 0(|i?|) time only: every element of R is moved from its old 
class Ci to its new class C^, the successor of Cj in the linked list. A flag in the data structure of Cj 
indicates wether its successor is or not. This allows to create C' if it does not already exist. C'^ 
is Ci n R while the remaining elements of Ci are \ R. All flags arc reset by a second scan of R, 
that also allows to remove empty classes from the partition hnked list (classes that were included 
in R). 

And at least, for a set Z and y ^ Z, Z can be partitioned according to the y-classes in Od.Z'j) 
time. If there are k y-classes, an array of k linked list is used and each element of Z is appended to 
the proper list. A stack of nonempty list allow to collect and reset them in 0(|2'|) time. 

Theorem 5. Algorithm 2 computes MHS{x) in 0{\V\^) time. 

Proof. For the correctness proof, one just has to check that the algorithm implements correctly 
Lemma 3. For the time complexity issues, notice that, for each pivot y, an element z is placed in 
Z only once. As partitioning Z into to y-classes, and then refining using all refining sets generated 
by y, takes 0(|.Z'|) time. Hence the algorithm takes 0(|F|^) time. □ 



7 



Algorithm 2: MMximid IIoiuogcuitxMiti Sets not coulaining x 

for every group G do 

for every class C of G do 

Compute the sot Z of elements of in G but not in C 
for every element y of C do 

Partition Z according to the y-classes 
|_ Add each partition set to the refining sets pool 

Set the group boundaries to the classes boundaries (from Pj-i to Pi) 
for each refining set R of the pool do 
|_ Remove R from the pool and then refine Pi using R 



6.3 Testing if a homogeneity relation is trivial 

A homogeneous relation i7 on F is trivial if 9'h contains only V and the singletons. 

Theorem 6. Let H be a homogeneous relation on V and S be a nonempty subset of V. One can 
test in 0(|yp) time if H is trivial. 

Proof. If |F| < 2 the answer is yes. Otherwise let x and y be two elements of F. In 0(|yp) time, 
Algorithm 2 outputs the maximal homogeneous sets not containing x. If one of them is nontrivial 
the answer is no. Otherwise all nontrivial homogeneous sets contain x. In 0(|yp) time, Algorithm 2 
outputs the maximal homogeneous sets not containing y. If one of them is nontrivial the answer 
is no. Otherwise all nontrivial homogeneous sets contain x and y. Then Algorithm 1 is used with 
S = {x, y}, in 0{\V\^) time. The answer is yes iff SHS{{x, y}) = V. □ 

6.4 Strong modules of a homogeneous relation 

Theorem 2 straightforwardly leads to an algorithm: 

Theorem 7. The strong homogeneous subsets of a homogeneous relation HonV can be computed 
in 0{\V\^) time. 

Proof. First compute MHS{x) for all x E V. All these sets together exactly form the family Z 
defined in Theorem 2. It can be done in 0(|yp) time using Algorithm 2 \V\ times. The size of 

this family (sum of the cardinals of every subsets) is 0(|y|^) since they form \ V\ partitions. Using 
Dahlhaus algorithm [7] the overlap components can be found in time linear in the size of the family, 
thus 0(|yp). According to Lemma 1 there are at most \V\ nontrivial overlap classes. For each class 
it is easy to compute its support, and in 0(|Fp) time easy to compute its atoms (each subset 
of the class is used as pivot, in a partition refinement of the support). And after all the 0(|yp) 
supports and atoms output must be sorted by inclusion order into the inclusion tree of the strong 
homogeneous sets (removing many duplicates), an easy task in 0(|yp) time. □ 

Notice that, if the homogeneous relation defines a weakly partitive family, then the quotient property 
applies and helps a lot. The algorithm scheme of [8], that can be implemented in 0(n+m) for graphs, 
could be implemented in 0(|y|^) time for an homogeneous relation. The approach is to compute 
MHS{x) then, using the quotient relation, to compute all strong homogeneous sets containing x. 
That gives the left branch of the decomposition tree. Then the algorithm is recursively launched. 
The amortised complexity analysis of Section 6.2 can be used: as a class is not splitted when the 
recursive process begins, the whole algorithm takes 0(|l^p). But our proof needs Axiom A^.- It 
holds for graphs, but not for directed graphs nor 2-structures, and seems very specific. 
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7 Applications 



Let us examine in the sequel some of the appUcations of this Homogeneity theory. 

7.1 Modular decomposition 

In a graph, the homogeneous relation H{x\yz) is true when x "sees" y and y in the same way. In 
undirected graphs, this means that either there are two edges (xy) and (xz), or no edge between 
X and the two other vertices. In directed graphs, this means that there are zero or two incoming 
arcs between x and the two other vertices, and zero or two out-coming arcs between x and the two 
other vertices. The homogeneous sets are then called modules. The notion of modules also extends 
to 2-structurcs [10]. A 2-structure is a complete edge-coloured graph and H{x\yz) is true when 
edges {xy) and (xz) have the same colour. 

Proposition 9. — The homogeneous relation of a undirected graph fulfills Al , A2, A3 and A4 

— The homogeneous relation of a directed graph fulfills A 2 and A3 

— The homogeneous relation of a 2-structure fulfills Al and A2 and A3 and A4 

The modules of a undirected graph and of a 2-structurc thus form a partitive family, while the 
modules of a directed graph just form a weakly partitive family. All know properties of modular 
decomposition [19] can be derived from this result. An 0{n^) modular decomposition algorithm 
can also be derived from Section 6.4 algorithm, but it is less efficient than the existing algorithms 
[14, 18,8, 16,6,4,13,9]. 

7.2 Other graph relations 

In a graph we can consider different homogeneous relations, for instance the relation "there exists 
a path from vertex x to vertex y avoiding the vertex s ", or a more general relation "there exists a 
path from x to y avoiding the neighbourhood of s ". It is easy to see that these two relations fulfill 
the basic axioms (symmetry, reflexivity and transitivity). In the first case, the strong homogeneous 
sets form a partition (into the 2- vertex-connected components, minus the articulation points). The 
second relation is related to decomposition into star cutsets. 

Another interesting relation is Dk{s\xy) if d{s,x) < k and d{s,y) < k, where d{x,y) denotes 
the distance between x and y. The case k = 1 corresponds to modular decomposition. It is worth 
investigating the general problem. 

7.3 Bimodular decomposition 

Let G = {B, W, E) be a bipartite graph where B contains the black vertices and W contains the 
white vertices. A bimodule is a subset of vertices M C {B U W) such that no black vertex b not 
in M distinguishes two white vertices of M (there must be either no or all possible edges between 
b and the white vertices of M) and conversely no white vertex not in M distinguishes two black 
vertices of M. In [11] is defined the bimodular decomposition of a bipartite graph. It is stated 
that, although the family of bimodules is not even partitive, then the strong bimodules are an 
optimal encoding of the family. Indeed, the inclusion tree of strong bimodules, plus some O(n^) 
pointers (easy to add, given the graph and the inclusion tree), are enough to store and output 
the (potentially exponential) family of bimodules, and to test if a set is a bimodule, and allow to 
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solve in polynomial time some NP-complete problems, when the degree of the nodes of the tree is 
bounded [15]. [11] give an O(mn^) time algorithm that compute the strong bimodules given the 
graph. But given the graph the homogeneous relation can be computed in O(n^) time by testing all 
refiectless triples, and then using the algorithm of Section 6.4 the strong bimodules can be output 
in O(n^) time, improving the previous time bound. 

8 Conclusion 

We hope that this homogeneity theory will have many other applications and will be useful to 
decompose automata [1] and boolean functions [2]. Obviously, the algorithmic framework presented 
here can be optimised in each particular application, as it can be done for modular decomposition 
[14, 18, 8, 16, 6, 4, 13, 9]. We think the homogeneity concept is a very general idea. 
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